Equilibrium structure and vibrational spectra of sila-adamantane 
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The recent synthesis of a four-fold silylated sila-adamantane molecule (C24H72Sii4, T^), [J. Fis- 
cher, J. Baumgartner, and C. Marschner, Science 310, (2005) 825] is the first attempt of making 
the silicon analogue of adamantane. It has adamantane-like Siio core, capped by methyl and sily- 
methyl ligands. We report its electronic structure, vibrational spectrum, and the infra-red and 
Raman spectra calculated within the density functional formalism using large polarized Gaussian 
basis sets. The properties of sila-adamantane are compared with exact silicon analogue (SiioHi4) of 
adamantane. Results show that replacing hydrogens in SiioHie by methyl and silymethyl ligands 
results in expansion of the Siio core and are responsible for large number of peaks in the Raman 
spectrum. The Si-C stretch at 664 cm~^ and methyl deformations frequencies compare well with 
recent measurements of vibrational frequencies of methylated silicon surface. 
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I. INTRODUCTION 

Silicon is of considerable interest in nanotechnology 
and in semiconductor industry. In bulk phase, silicon 
has diamond lattice like carbon. The diamond structure 
can be built by translating the smallest unit of the car- 
bon lattice called adamantane. The adamantane has 10 
carbon atoms that are capped by hydrogen atoms that 
passivate the dangling bonds of carbon atoms in the sp^ 
bonding. The adamantane and a few low-order diamon- 
doids have been synthesized [1]. Recently, isolation of 
several new lower- and medium-order diamondoids from 
petroleum oil was reported by the Molecular Diamond 
Technologies Group at Chevron- Texaco [2|. This has led 
to a renewal of interest in the chemistry and physics of 
the diamondoids [3]. 

The search for silicon analogue of adamantane, called 
sila-adamantane SiioHig, has so far been unsuccessful. 
There has been however some progress. A recent report 
by Ficher, Baumgartner and Marschner (FBM) of syn- 
thesis of silylated sila-adamantane is a breakthrough in 



the synthesis of sila-adamantane^ . Instead of hydrogen 
atoms as capping elements as in the exact silicon ana- 
logue of adamantane, silicon atoms in sila-adamantane 
are capped by methyl and trimethylsilyl groups (Cf. Fig. 
[1} . It was characterized using the nuclear magnetic reso- 
nance (NMR) spectroscopy. The ultra-violet and visible 
spectra were also measured which showed strong absorp- 
tion at 222 nm (5.58 eV). Followed by this report, Pic- 
cheri reported electronic structure of C24Sii4H72 within 
Perdew-Burkc-Ernzerhof density functional model. He 
reported structural parameters calculated using the 3- 
21G* basis sets and density fitting methodology [5|- 

In the present article, we report the electronic struc- 
ture of sila-adamantane examined within the density 
functional theory using large polarized Gaussian basis 
sets. Furthermore, at the same level of theory we cal- 
culate vibrational frequencies, the infra-red and Raman 
spectra. We also calculate these properties for exact 
silicon analogue of adamantane (SiioHig), and compare 
them with those of sila-adamantane C24Sii4ll72 to ex- 
amine the effect of methyl and silymethyl ligands. 




FIG. 1: (Color online) The optimized structures of 
sila-adamantane SiioHig and silylated sila-adamantane 

C24Sil4H72. 








II. FORMALISM 

Our calculations are performed within the Kohn-Sham 
formulation [6| of the density functional theory 7| us- 
ing NRLMOL code'!l|. The molecular orbitals are ex- 
pressed as linear combination of Gaussian functions. 
The exchange-correlation effects are described within the 
generalized gradient approximation using Perdew-Burke- 
Ernzerhof[8J parameterization. For a reliable prediction 
of properties, we use large polarized basis sets that are 
specially optimized for the PBE functional[10|. Thus the 
basis for Si contains 6 s-, 5 p-, and 3 d-type Gaussians 
each contracted from 16 primitive functions. For C, 12 
primitive functions were used to contract 5 s-, 4p-, and 3 
d-type basis functions. The basis for H consist of 4 s-, 3 
p-, and 1 d- type basis functions. In total 2614 basis func- 
tions were used in calculation of C24Sii4H72. Additional 
diffuse functions were used for the calculations of polar- 
izabilities. The exchange-correlation contribution to ma- 
trix elements and total energy are obtained numerically 
using efficient numerical grids employed in NRLMOL [9j. 
The vibrational frequencies were obtained by diagonaliz- 
ing the dynamical matrix built from the forces computed 
at perturbed geometries. Finally, the infra-red and Ra- 
man spectra were computed from the dipole moment and 
polarizability derivatives, respectively. 



FIG. 2: (Color online) The optimized structure of C24Sii4ll72 
without hydrogens. The Si atoms that are part of the sily- 
methyl ligands are shown in dark blue for clarity. The atomic 
labels are identical to those in experimental structure from 
Ref. y. 

III. RESULTS 

The optimized structure of SiioHig and C24Sii4H72 
are shown in Fig. 1. The C24Sii4H72 is also shown with- 
out hydrogens in Fig. 2. To distinguish the Si atoms that 
belong to silymethyl-ligand, these four atoms are shown 
in dark blue. The SiioHig core in C24Sii4H72 is clearly 
seen. The atoms have labels exactly identical to those 
given in the supplementary information of the experi- 
mental study in Ref. |j. Both structures were optimized 
within the T^ symmetry. The bond distance between two 
silicon atoms (Si(l)-Si(2)) that belong to the tetrahedral 
adamantane core (Si(lO)) is 2.39 A, in good agreement 
with experimental value of 2.36 A. In case of SiioHig, 
this bond length is calculated to be 2.36 A. The Si-SiCHa 
bond length (Si(2)-Si(5) distance) is also 2.39 A while the 
experimental value is 2.35 A. The bond angle between 
the Si-Si bonds from the Siio core (Si(2)-Si(l)-Si(2A)) in 
C24Sii4H72 is 110° and is in excellent agreement with the 
experimental value of 110.3°. Similarly, very good agree- 



ment is observed for the dihedral angle between three 
bonds from the Siio core ( Si(2A)-Si(l)-Si(2)-Si(3)). The 
predicted angle is 59.8° while experimental value is 59.9°. 
These bond angles are essentially similar in SiioHig (an- 
gles 110.1° and 59. 7°). As the bond angles and dihedral 
angles in SiioHig are essentially similar to those in the 
Siio core in C24Sii4H72, it appears that replacing hy- 
drogens by methyl and trimethylsilyl ligands leads to an 
expansion of the Siio core. 

IV. ELECTRONIC STRUCTURE PROPERTIES 

The fully occupied highest occupied molecular or- 
bital (HOMO) of C24Sii4H72 is three fold degenerate 
and has T2 symmetry while its LUMO has Ai sym- 
metry. The HOMO-LUMO (HL) gap is 4.16 eV. For 
the SiioHie, the HL gap is 4.82 eV. The HL gaps re- 
ported earlier with a smaller 3-21G* orbital basis and 
the density fitting methodology are in the range 4.38- 
4.53 eV depending on the choice of exchange-correlation 
parameterization [5|. We also calculated the HOMO- 
LUMO gap for the SiioHig using the time dependent 
DFT (TDDFT) and using the Perdew-Wang GGA as 



implemented in ABINIT ll| . The calculated gap using 
TDDFT of 4.6 eV is only marginally larger than the 
PBE-DFT gap. The optical gaps of SiioHig have been 
studied using more complex quantum Monte Carlo and 
GW approximations [1^. These methods give value of 9.2 
eV for the optical gaps. Our attempts to compute the 
TDDFT optical gap for C24Sii4H72 were not successful 
due to large memory requirement. 



V. VIBRATIONAL SPECTRA 

A. SiioHig 

The SiioHig cage has T^ point group symmetry. It has 
in total 72 vibrational modes which could be classified as 
5Ai -I- A2 + 6 E + 7Ti + 11 T2. Of the total 11 T2 
modes that are infra-red active, only four have apprecia- 
ble intensity. The calculated IR spectra is shown in Fig. 
H The first absorption is at 716 cm"! (7.61 D/A^) due 
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FIG. 3: (Color online) The infra-red spectra of SiioHie and 

C24Sil4H72. 
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FIG. 4: (Color online) The Raman spectra of SiioHig. 
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FIG. 5: (Color online) The Raman spectra of C24Sii4H72. 
It is split into low frequency region 50-1000 cm^^ and high 
frequency 2800-3100 cm~^ (CH stretches) region. 



to the SiHa wagging. The mode at 885 cm-i (3.85 D/A^) 
corresponds to the bending (scissor) mode of SiH2. The 
peaks for the two SiH stretching modes are at higher fre- 
quencies of 2130 cm~^ and 2141 cm~^ are coalesced in 
the figure. The quantities in bracket are intensities. 

The Raman spectra of SiioHig is shown in Fig. U) 
There are 5Ai, 6 E, and 11 T2 Raman active modes. 
The fully symmetric Ai mode at 2137 cm^^ is the sym- 
metric SiH stretch and is the most intense (intensity 130 
A'*/amu). The absorption due to fully symmetric stretch 
of Si-Si bonds (breathing of silicon cage) occurs at lower 
(due to larger mass) frequency of 299 cm~^ . The second 
most intense absorption is seen at lower wave number of 
102 cm~^ with intensity of 68 A^'/amu. This is three fold 
degenerate Ti mode due to rocking motion of SiH2. The 
Raman active doubly degenerate mode at 103 cm^^ is 
due to the deformation of silicon cage. The absorptions 
at larger frequencies ( > 600 cm^^ ) involve SiH2 defor- 
mations and SiH stretches. Less intense absorption are 
620 cm~i (SiH2 twist, intensity 9 A**/amu), 688 cm~i 
(SiH2 and SiH rocking, intensity 2 A'^/amu), and 878 
cm~^ (SiH2 scissor, 10 A^/amu ). 



B. C24Sil4H72 

We computed the vibrational frequencies, within the 
harmonic approximation, with different possible orienta- 
tions of methyl groups that are consistent with the Td 
symmetry. In all cases we find 9 imaginary frequen- 
cies for C24Sii4H72 due to the rotation of CH3 group. 
These rotations seen in the gas phase will be frozen in 
the crystalline state. The total 324 vibrational modes in 
the C24Sii4H72 can be classified as 17Ai -I- 10 A2 + 27 
E + 37Ti + 44 T2. There are a total of 44 IR active 
modes. The calculated infra-red absorption spectrum is 
displayed in Fig. |3l The two conspicuous peaks at lower 
frequencies of 795 cm^^ (8 D/A^) and 834 cm"^ ( 9 
D/A^) are due to the SiCHa rocking/wagging motions. 
These modes are similar to the SiH2 wagging and rock- 
ing seen in SiioHigThe weak absorption at 1232 cm~^ 
( 1.4 D/A^) is the Si-CHa umbrella mode. The three 



peaks at higher frequencies in the frequency range 2956- 
3032 cm~^ with intensities of 1.3-1.4 D/A^ are due to 
the symmetric and asymmetric C-H stretches. A number 
of vibrational modes results in weak absorption, with in- 
tensities less than 1 D/A^. These are Si-C asymmetric 
stretch at 664 cm~^ (0.6 D/A^), Si-Si-Si deformation 
at 631 cm^i (0.7 D/A^), Si-C symmetric stretch at 603 
cm-i (0.6 D/A2), CH3 deformation at 1417 cm^^ (0.7 
D/A2), and C-H symmetric stretch at 3060 (0.8 D/A^). 
The compilation of the frequencies is given in Table HI 

TABLE I: The IR active frequenters for the C24Sii4H72. The 
stretch of absorptions is given in bracket with following ab- 
breviations: w: weak, vs: very strong, m: medium. The 
experimental numbers are for methylated silicon(lll) sur- 
face obtained by high-resolution electron energy loss spec- 
troscopy (HREELS) [131 and transmission infra-red (TIR) [141 ] 
spectroscopy. 



Freq. ( cm 


"^ ) HREELS 


TIR Mode description 


664 (w) 


678 


Si-C asymmetric stretch 


603 (w) 




Si-C symmetric stretch 


631 (w) 




Si-Si-Si deformation 


796(vs) 




SiCHa bending 


834(vs) 




CII3 rocking 


1233(m) 




1257 CH3 Umbrella 


1417 (w) 


1424 


2856 CH3 deformation 


2956 (m) 




2909 C-H symmetric stretch 


2969 (m) 


2916 


2909 C-H symmetric stretch 


3032 (m) 


2989 


2928 C-H asymmetric stretch 



The Raman spectrum of C24Sii4H72 is shown in Fig. 
[5] The spectrum can be split into two parts: the ab- 
sorptions in the low-intermediate frequency range (60- 
200 cm^^ ) and intense absorptions in the frequency 
range 2900-3100 cm~^ . The doubly degenerate absorp- 
tion at 59 cm"^ (77 A*/amu) involves rocking like mo- 
tion of Si-Si(CH3)3. This mode is depicted in Fig. [51 
The absorption at 84 cm~^ (63 A^/amu) is due to a 
doubly degenerate mode that results in deformation of 
Si-Si-Si bonds in the Siio core similar to the mode at 103 
cm~^ seen in SiioHigThe triplet degenerate mode also 
due to the distortion of Siio core at 85 cm~^ is merged 
with the peak at 84 cm^^ . These two modes are shown 



in Fig. [T] The weak absorptions at 108, 155 cm~^ are 
due to methyl rotations. The most intense peak in the 
low frequency region at 130 cm"-'^ (85 A'*/amu) is due 
to non-degenerate Ai mode that resuhs in breathing mo- 
tion of Si atoms to which Si(CH3)3 ligands are attached. 
The breathing motions of rest of Si atoms in Siio core (to 
which two CH3 groups are attached) leads to an absorp- 
tion at 177 cm~^ ( A'^/amu). 

The intense absorptions at higher frequencies are due 
to CH stretches. The first peak at 2956 cm~^ (160 
A'*/amu) is due to Ai CH stretch where CH3 belongs 
to the Si(CH3)3 group. The Ai CH stretch in the two 
methyls attached to silicon atoms leads to absorption at 
slightly higher frequency 2969 cm~^ . The three- fold de- 
generate asymmetric CH stretches absorb at 3032-3033 
cm^^ ( 102 A^/amu) and is also infra-red active. 

The most obvious difference between the Raman spec- 
tra of SiioHig and C24Sii4H72 is the presence of the 
intense peak at 2130 cm^^ for the SiioHig involving 
SiH bonds and the high frequency modes involving CH 
bonds in the methylated C24Sii4H72 cluster. It may be 
mentioned here that there are several very low inten- 
sity modes between 1000-1500 cm~^ for the C24Sii4H72. 
However, the intensities of these peaks are comparatively 
very small and hence may not be used for identification 
purposes. 

The passivation of dangling bonds at the silicon sur- 
faces by hydrogenation, chlorination and alklation (ethy- 
lation, methylation etc) has been a subject of sev- 
eral studies due to technological importance of silicon. 
Two groups have reported experimental studies of vi- 
brational modes of methylated silicon surface have been 
reported[lJ|. Although the sila-adamantane core studied 
here is a very small unit of silicon crystal, it is inter- 



esting to compare the frequencies due to Si-C stretching 
and methyl deformations. Yamada and coworkers 13 as- 



sign the vibration of 678 cm~^ to Si-C which compare 

well with 664 cm^^ Si-C asymmetric stretch. The CH3 

stretch and deformation freuenqcies compare well (Cf. HI. 

VI. CONCLUSION 

The equilibrium structure of recently synthesized 
of silylated sila-adamantane C24Sii4H72 was obtained 
within the density functional theory, using large polarized 
Gaussian basis sets and the semi-local PBE approxima- 
tion to the exchange-correlation functional. The calcu- 
lated geometrical parameters are in very god agreement 
with experimental values. The Siio core in the optimized 
C24Sii4H72 is dilated with respected to the silicon ana- 
logue of adamantane SiioHig. The calculation of vibra- 
tional frequencies indicate that some of the methyl lig- 
ands are freely rotating in the gas phase. The infra-red 
and Raman spectra of C24Sii4H72 is reported for the first 
time. The comparison of spectra for the C24Sii4H72 with 
that of SiioHie do not show any similarity since in both 
the structures only the Si core is similar. The SiioHig 
shows intense peaks corresponding to stretching of the 
Si-H bonds which are missing in the C24Sii4H72 due to 
absence of SiH bonds. On the other hand, the presence 
of methyls are evident in the low frequency peaks for the 
C24Sii4H72 and also in very high frequency peak arising 
from the CH bonds. 
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